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FREQUENCY RESPONSE AND TRANSFER FUNCTIONS OF A NUCLEAR ROCKET 
ENGINE SYSTEM OBTAINED FROM ANALOG COMPUTER SIMULATION 
by Clint E. Hart and Dale J. Arpasi 
Lewis Research Center 

SUMMARY 

The dynamic response characteristics of the basic variables of a NERVA-type nuclear 
rocket engine system are presented. Frequency-response data, at four selected opera- 
ting levels, were obtained from a detailed analog computer simulation. Frequency re- 
sponses of thrust chamber temperature, chamber pressure, and reactor power to perti- 
nent control or system variables are shown. 

Transfer functions that approximate the frequency responses are presented. The ef- 
fect of reactor period on the reactor -power transfer function was investigated, and trans- 
fer functions are presented for 2. 0 and 0. 5 seconds, as well as for infinite periods. These 
results provide essential information for control system analysis and design. 

The analytical models and assumptions used to develop the equations that represent 
the dynamics of nuclear rocket engine components are discussed. A set of specific equa- 
tions that can be used in an analog computer simulation is presented. 


INTRODUCTION 

Extensive experimental testing is being conducted in the NERVA program to provide 
data for the design of a nuclear rocket engine system. Although the type of engine system 
under consideration was found to operate stably without closed-loop control, transient 
operational requirements necessitate the use of a closed-loop control system. These 
requirements include rapid startup (possibly 60 to 90 sec from flow initiation to full 
thrust), precise control of thrust and specific impulse, and controlled shutdown and re- 
start . 

For the design of such a control system, either transfer functions or frequency re- 
sponses of some basic engine system variables to other system or control variables are 
required. During the experimental tests conducted thus far, two methods were used to 



obtain these frequency responses or transfer functions. These methods were (1) direct 
measurement of the frequency responses to sinusoidal input signals and (2) measurement 
of responses to random noise input signals and use of cross-correlation data processing 
techniques (ref. 1). 

Neither method yielded sufficient frequency-response data or transfer-function infor- 
mation to complete a comprehensive control-system design. The main reasons for the 
scarcity of this information are the difficulty of obtaining measured signals with good 
signal-to-noise ratios and the limited amount of running time allotted to this type of test. 

System transfer functions can be analytically derived from the set of algebraic and 
nonlinear differential equations that describe the system behavior (refs. 2 and 3). How- 
ever, this is a formidable task. Many simplifying assumptions are required to derive 
reasonable transfer functions . The validity of these assumptions has not been thoroughly 
investigated. 

A fourth method is to simulate the nuclear rocket engine system by programing an 
analog computer to solve a set of system equations and to obtain the required frequency- 
response data from this simulation. 

The purposes of this study were (1) to obtain frequency responses of basic engine 
system variables to other system or control variables at several steady- state operating 
levels and (2) to determine the approximate transfer functions that represent these fre- 
quency responses. 

A set of equations was developed that represents the principal dynamic characteris- 
tics of a nuclear rocket engine system over a range from approximately 20 to 120 percent 
of the design conditions of thrust chamber pressure and temperature. A nonlinear, wide- 
range analog computer simulation was programed from this set of equations; frequency 
responses of reactor power, thrust chamber pressure and temperature were obtained at 
four selected operating levels. 

Approximate transfer functions that represent these frequency responses were de- 
termined. In addition, reactor-power transfer functions that are dependent on reactor 
period were developed. These transfer functions are essential for control system de- 
sign by root-locus techniques. 


SYSTEM CONFIGURATION 

The nuclear rocket engine system considered in this study is shown schematically in 
figure 1 . The pump receives the liquid hydrogen, that is stored at low pressure in the 
propellant tank and delivers it at high pressure through a discharge line to the inlet of the 
nozzle coolant tubes. The hydrogen vaporizes in the nozzle coolant tubes. Then the hy- 
drogen vapor flows through the reflector, shield, and reactor core passages and is ex- 
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Figure 1. - Schematic drawing of nuclear rocket engine. 










pelled at high velocity through the nozzle to provide thrust. Heat is transferred, prin- 
cipally by convection, to the hydrogen as it flows through the system. 

The hot -bleed cycle provides the turbine hot-gas source. Gas of the proper temper- 
ature to drive the turbine is obtained by controlled mixing of hot gas tapped from the reac- 
tor outlet chamber and cold gas from the reactor inlet plenum. Turbine weight flow, and 
hence torque, is controlled by a valve in the bleed line. 

The reactor is a solid-core thermal type with highly enriched uranium fuel and a 
graphite moderator. The reflector is a beryllium sleeve surrounding the core. Cylin- 
drical control drums, with neutron poison plates covering part of their surface, are ro- 
tated to control reactor power. These control drums are located in the reflector. 

SYSTEM EQUATIONS AND ASSUMPTIONS 

The equations that describe the dynamic behavior and interrelation of the engine sys- 
tem variables consist of a set of ordinary and partial differential equations and some 
algebraic equations. To develop a set of equations for analog computer simulation of a 
system, the geometric complexities of the system must be reduced and a simpler, ana- 
lytical model formulated. Furthermore, the partial differential equations must be re- 
duced to ordinary differential equations since the analog computer normally performs 
operations with respect to only one independent variable. Finite -difference techniques 
were used to reduce the equations to time -dependent differential equations. 

In the following paragraphs of this section, the analytical models and assumptions for 
the various components comprising a nuclear rocket engine system are discussed and the 
general equations that evolve from these models are given. The specific equations used 
in the analog computer simulation are given in appendix B. These equations are based on 
the design data for a first-generation NERVA-type nuclear rocket engine that develops 
approximately 56 000 pounds of thrust. The station designations used in the specific 
equations are shown in figure 1 . 


Reactor Core 

The analytical model assumed for developing the reactor -core heat-transfer equa- 
tions is a cylinder with a single flow passage divided axially into three equal-length sec- 
tions. The following discussion applies to one section. 

Taking a heat balance on the material and assuming the average wall temperature to 
be equal to the average material temperature yielded the following equation: 
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( 1 ) 
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dt 


= Q m - hA c 
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[ T T in + T out 
\ m ~ o 


(All symbols are defined in appendix A.) The first term on the right side represents the 
thermal power generated in the section, and the second term represents the heat trans- 
ferred by convection to the propellant. The heat content of the core section was obtained 
by integrating equation (1) . _ 

The average material temperature T m is considered to be a function of the heat 
content q m< Since the specific heat of the core material varies considerably with tem- 
perature, the heat content is defined as the mass times the integrated specific heat of the 
material, and is in equation form 


q 


m 


= p V 
'm m 



m 



( 2 ) 


The heat -transfer coefficient h in equation (1) was determined by using the following 
correlation equation from reference 4 


0.025 Re 0 ’ 8 Pr 0 ’ 4 



(3) 


with k, Re, and Pr evaluated at bulk conditions. The L/D correction was negligible 
for the analytical models considered in this study. 

The following heat balance equation was obtained for the propellant contained within 
the flow passage: 


dt 


= hAj T 


s\ m 


T in + T out\ 


W < H out - H in> 


(4) 


In the system under study, the thermal time constant of the propellant is of the order of 
1 millisecond, while the thermal time constant of the core material is of the order of 
1 second. Thus, in the frequency range of interest (up to 100 rad/sec), the thermal dy- 
namics in the propellant can be neglected. Equation (4) then can be reduced to 


hA 


T T in + T out] 


s\ m 


= W < H out * H in> 


(5) 
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The enthalpy of the propellant at the core inlet is considered to be a function of pres- 
sure and temperature at temperatures below 240° R, but only a function of temperature at 
temperatures at or above 240° R. The enthalpy of the propellant at other locations in the 
core is considered to be only a function of temperature. 

The momentum or pressure -drop equation used in this study is 


P in 



v in + 


4fL/l 


D\p 


+ L dW 
144gA dt 


( 6 ) 


Over the range of conditions encountered in the flow passages, the term 1/p can be 
approximated by (v- n + v ou p/4. Inlet and outlet specific volumes were computed from 
pressures and temperatures by use of the state equation for a perfect gas. Additional 
pressure losses at the inlet and exit of the flow passage were also considered. 

When the momentum equation (eq. (6)) was applied to the reactor core, the dynamic 
contribution of the last term was assumed to be negligible, and the equation was solved 
for weight flow rate. As a compromise between accuracy and the analog components re- 
quired, only one axial section was assumed for the reactor core when equation (6) was 
used. The friction factor was assumed to be constant in the reactor core. 


Reactor Kinetics 

The reactor-kinetics equations were used in this study in the standard linear form as 
follows: 


dn= 5K n -An + f A.C. 
dt l* l* fr[ 1 1 


( 7 ) 


dC. fi. 

= — n - A-Cj 

dt l* 11 


( 8 ) 


The values of the decay constants and the relative abundance of the six delayed-neutron 
groups were obtained from reference 5. 


Reactivity 


The total reactivity increases or decreases reactor power according to the reactor- 
kinetics equations. Total reactivity consists of the inherent reactivity feedback plus the 


6 



reactivity due to the control drums. The inherent reactivity feedback is a rather compli- 
cated phenomenon, which, for the purposes of this study, was assumed to consist of the 
reactivity related to material temperature changes in the core and to hydrogen density 
changes in the core flow passages. 

The inherent reactivity feedback is given by the equation 

6K l' F < ? mc’PH> < 9 > 


The average material temperature in the core T mc is approximated by the linear aver- 
age of the material temperatures in the three axial sections. The density of the hydrogen 
in the core changes over a wide range from inlet to outlet and is not linear with core 
length. Thus, the average hydrogen density Pjj cannot be determined by simple aver- 
aging techniques. For the range of variables encountered in this study, the following 
equation is a good approximation of the average hydrogen density in the core: 


P H = 


144P in 


R(T 


out 


T. ) 
in' 


In 


out 

T. 

in 


( 10 ) 


Pump 


The pump inlet pressure and temperature were assumed to be constant and equal to 
tank conditions. The pump pressure rise and torque were determined from functional re- 
lations among pressure rise, torque, weight flow rate, and speed. Typical pump maps, 
which were obtained from experimental data by plotting pressure rise against weight flow 
rate for constant speeds and constant efficiencies, were converted to single-curve func- 
tional relations defined by the following equations : 



Density was assumed to be constant throughout the pump. The temperature rise be- 
tween the pump inlet and outlet was assumed to be proportional to weight flow for this 
study. 
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Pump Discharge Line 


The momentum equation (eq. (6)) was solved for dW/dt and integrated to obtain 
pump weight flow rate. The specific volume and the friction factor were assumed to be 
constant, and the temperature rise in the discharge line was assumed to be zero. 


Nozzle 

The analytical model assumed for the nozzle coolant passages was a single tube with 
one axial section. Although the tube geometry was quite complex, equivalent constant- 
area tube properties were used in the pressure-drop and heat -transfer equations. The 
momentum equation (eq. (6)) was solved for the inlet pressure; the dW/dt term was 
neglected. Inlet specific volume and friction factor were assumed to be constant. Outlet 
specific volume was determined from outlet pressure and temperature by using the state 
equation. 

Nuclear heat generation in the nozzle wall material and thermal radiation to the noz- 
zle walls from the hot core were assumed to be negligible. Thus, a heat-balance equation 
for the nozzle wall material can be written as 

~ = WV ' ? m> - \c A s - T ‘ n * 2 T ° Ut ) <»■ 3 > 

The first term on the right side represents the heat transferred from the hot gas to the 
nozzle wall material, and the second term represents the heat transferred from the ma- 
terial to the hydrogen flowing through the coolant tubes. In this study the nozzle heat- 
transfer coefficients are assumed to vary only with weight flow to the 0. 8 power. Equa- 
tion (13) was integrated to obtain the heat content of the wall material. As in the core, 
the material temperature was considered to be a function of the heat content. 

The nozzle coolant outlet enthalpy and temperature were obtained by solving equa- 
tion (5) for outlet enthalpy and considering outlet temperature as a function of enthalpy and 
pressure. 


Reflector 

Although the reflector contains different materials and also flow passages of various 
sizes and shapes, it is not feasible in this system -dynamics study to solve equations rep- 
resenting the heat transfer and fluid flow of all the different types of flow passages in 
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parallel. Instead, a single cylinder and flow passage with one axial section was assumed 
for the reflector model. The size of the flow passage was determined from total flow area 
and heat-transfer surface area. The physical and thermal properties of the cylinder ma- 
terial were based on weighted averages of the properties of aluminum, beryllium, and 
graphite. 

The thermal dynamics in the reflector are represented by equations (1) to (3) and (5) 
with the proper coefficients. Equation (5) was solved for coolant outlet temperature. The 
enthalpy of the hydrogen at the reflector inlet was considered to be a function of both tem- 
perature and pressure because, in this study, the temperature at this point was always be- 
low 240° R. The reflector pressure drop is represented by the momentum equation 
(eq. (6)). However, this equation was solved for the inlet pressure and the dW/dt term 
was neglected. 


Reactor-Core Inlet Plenum and Shield 

The density of the hydrogen in the core inlet plenum was obtained by integrating the 
following continuity equation: 


ydP=Vw (14) 

dt ^ 

The pressure in the plenum was computed by using the state equation and assuming uni- 
form density and temperature throughout the plenum . The heat transfer and pressure 
drop in the shield were considered to be negligible for this study . 


Thrust Chamber and Nozzle 

The continuity equation (eq. (14)) and the state equation were used to determine the 
density and pressure of the hydrogen in the thrust chamber. For this study, the thrust 
nozzle was assumed to be choked at all times and nozzle weight flow rate was computed 
from the equation 


W = K 


P 

rpO. 5 


(15) 
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Bleed System and Turbine Control Valve 


A portion of the turbine bleed line was considered to be a mixing volume. The den- 
sity and pressure in the mixing volume were computed by using the continuity equation 
(eq. (14)) and the state equation. The hot bleed gas was considered to flow through a 
fixed orifice and the cold bleed gas through a valve in the cold bleed line. 

The cold-bleed, hot-bleed, and valve weight flow rates were computed by using the 
following compressible -flow equation: 


W = A 


p 

2g| 

f r \ 

rpO. 5 

_R ^ 

V - l/_ 


-|0.5 


(PR) 1 /^ - (PR) (r-1) / r ] 


0.5 


(16) 


In this equation, P and T are values upstream of the orifice or valve, and PR is the 
ratio of downstream to upstream pressures. The turbine control valve is a butterfly -type 
valve. The effective flow area is a function of the shaft angular position of the valve. 

Since the cold-bleed, hot-bleed, and valve weight flow rates were at different temper- 
atures, a form of energy -balance equation was used to compute the heat content of the hy- 
drogen in the mixing volume. The enthalpy of the hydrogen in the mixing volume was 
computed by integrating the following equation: 

PV ^ = w hb H hb + *cb H cb - w v H dV 

Then, from an enthalpy -temperature relation, the average temperature in the mixing vol- 
ume was obtained. To apply this equation in this study, the dynamics associated with 
mass storage were assumed to be negligible. 


Turbine and Turbopump Dynamics 


The section of line from the turbine control valve to the turbine inlet was considered 
the turbine inlet plenum. The turbine inlet temperature was assumed to be equal to the 
temperature in the bleed mixing volume. The continuity equation (eq. (14)) and the state 
equation were used to determine the density and pressure in the turbine inlet plenum. 

The turbine weight flow rate was computed from the following equation: 



( 18 ) 
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For pressure ratios less than 0.3, the functional relation is essentially constant. The 
turbine torque was computed by using the following equation: 



L - 1 J 


The pressure ratio used in equations (18) and (19) is the ratio of turbine-exit static pres- 
sure to turbine -inlet total pressure. 

An outlet plenum and exhaust nozzle were assumed to be downstream of the turbine. 
The density in the outlet plenum was determined by integrating the continuity equation 
(eq. (14)). The weight flow rate through the exhaust nozzle was computed by using the 
compressible -flow equation (eq. (16)) with appropriate constants. 

The enthalpy drop across the turbine was determined from the power output of the 
turbine (i.e. , the product of torque and speed). Then, if specific heat is assumed to be 
constant, the following equation for turbine exit temperature can be derived: 


T 


te 


M t N7T 
30W.C J 

L p 


( 20 ) 


In this study, the static temperature computed by this equation was assumed to be equal to 
the total temperature. 

The total pressure in the turbine outlet plenum was computed by using the state equa- 

2 

tion. The static pressure can be determined by subtracting pu /2g from the total pres- 
sure. The resulting equation for turbine-exit static pressure is 


Pte * P te 


RT te W t 

2 gP te A 2 


( 21 ) 


The turbopump speed was obtained by integrating the following torque -balance equa- 
tion: 


— = 30 ( M _ MJ (22) 

dt jrl 1 p 


11 



COMPUTER SIMULATION OF SYSTEM 
Analog Computer Simulation 


The system equations listed in appendix B were programed and set up on two 100- 
amplifier analog computer consoles. The complement of computer components required 
consisted of 190 amplifiers, 70 multipliers, and 20 diode -function generators. The time 
scale was 1 to 1 (i. e. , real time). Amplitude scaling permitted the simulation to cover 
the following ranges of the primary variables: temperature, 50° to 5000° R; pressure, 
10 to 1000 pounds per square inch absolute; weight flow rate, 5 to 100 pounds per second; 
and reactor power, 20 to 2000 megawatts. 


Digital Computer Solution of Steady-State Equations 

Amplitude scaling a large analog simulation to operate over a wide range of the sys- 
tem variables is a difficult problem. Optimum scaling at one operating point may cause 
scaling problems at other points. Knowledge of the steady -state values of the system 
variables at many operating points is required to scale the analog simulation properly for 
maximum accuracy over a wide operating range. Therefore, the steady-state system 
equations were programed for digital computer solution. 

The steady-state values obtained from the digital computer program were used not 
only to assist in scaling the analog simulation but also to check the validity of the analog 
setup and the overall accuracy of the analog simulation. 

One way to present the results of the steady-state calculations is shown in figure 2. 


8 . 

£ 


2 

E 


5000 


4000 


3000 


2000 


1000 ' 



200 300 

Chamber pressure, psia 

Figure 2. - Map of steady-state operating conditions with constraints. 
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Lines of constant control-drum reactivity and constant turbine -control-valve position are 
plotted as functions of chamber temperature and pressure. The boundaries or constraints 
on the allowable operating region are imposed by system component limitations or opera- 
tional requirements. The sensitivities of engine system variables to changes in the input 
variables can be estimated from this figure, but more complete information regarding 
system dynamics is presented and discussed in the following sections. 


FREQUENCY RESPONSE AND TRANSFER FUNCTIONS 

An analog computer simulation such as the one described herein can be used to obtain 
the frequency response of the engine system variables to changes in the input variables. 
This information is needed by the control-system designer to select feedback control-loop 
configurations and controller compensations properly, and to specify dynamic perfor- 
mance requirements of the control -loop hardware. 

One control concept that has evolved from various studies is, basically, to control 
reactor outlet -chamber temperature by manipulating the control drums and reactor outlet- 
chamber pressure by manipulating a turbine control valve. Figure 3 is a block diagram 
of this control concept. In addition to the primary control loops, a reactor -power control 
loop is used as an inner loop to the chamber temperature loop. In this control concept, 
the primary frequency -response data pertinent to control-system design include 



Figure 3. - Block diagram of control system for nuclear rocket engine. 
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(1) Response of chamber temperature to reactor power, ATg/AQ 

(2) Response of chamber pressure to turbine -control- valve position, APg/A0 v 

(3) Response of reactor power to control-drum reactivity, An/ASKj^ 

The control system may be required to provide stable operation over a wide range of 
the system variables. Therefore, frequency -response data were obtained from the ana- 
log computer simulation at design level and three other operating levels. The four oper- 
ating levels are shown in figure 2. The steady-state values of reactor power, core 
weight flow, chamber pressure and temperature are tabulated in table I. 

TABLE I. - STEADY-STATE VALUES OF SYSTEM VARIABLES 


Steady-state 

operating 

level 

Thermal power 
generated (ocn), 

Q> 

Btu/sec 

Reactor 

outlet 

chamber 

temperature, 

T S’ 

°R 

Reactor outlet 
chamber total 
pressure, 

P 8’ 

psi 

Reactor 
weight 
flow rate, 

W c> 

lb/sec 

I 

1.067X10 6 

4090 

550 

72.15 

n 

3.252X10 5 

2000 

253 

48.03 

in 

2.122 

2000 

167 

31.43 

IV 

4.780 

4090 

250 

32.41 


Frequency-response data can be used directly to determine graphically, by means of 
Bode plots and Nichols charts, the open-loop and closed-loop frequency response of the 
system with various feedback control-loop configurations and controller compensations. 
Another graphical analysis technique that can be used is the root-locus technique. This 
technique generally is faster and provides more insight into the effects of the controller 
on the transient response of a system. However, to apply the root-locus technique, it is 
necessary to have the system dynamic response in transfer -function form. 


Equipment and Procedure 

A frequency -response analyzer was used to obtain the frequency-response data. This 
equipment consists of a variable -frequency sinusoidal oscillator and a signal analyzer. 

The analyzer calculates the amplitude ratio and phase shift between an input and output 
signal of the component or system being tested or simulated. 

With the analog computer simulation operating at a steady-state level, a sinusoidal 
signal from the oscillator was added to one of the input variables. The resulting response 
of the selected output variable was fed to the signal analyzer, which calculated and dis- 
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played an amplitude -ratio and a phase-shift reading. The amplitude -ratio reading was 
corrected to account for amplitude scaling of the analog computer signals. This proce- 
dure was repeated for other input signal frequencies over the range of interest and also 
for various combinations of input and output variables. 

An Electro Scientific Industries Algebraic Computer (ESIAC) was used to determine 
transfer functions that approximate the frequency -response data. This special type of 
analog computer and some of its applications are described in references 6 and 7. For 
this application, the poles and zeros of the trial transfer function were adjusted until the 



(b) Phase shift. 

Figure 4. - Frequency response of reactor-outlet-chamber temperature to reactor power, obtained 
from analog computer simulation. Turbine-control-valve position held constant at each steady- 
state operating level. 
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frequency response calculated by the ESIAC computer matched the analog frequency- 
response data to the desired accuracy. An agreement of ±5° in phase shift and ±10 per- 
cent in amplitude was considered satisfactory. 


Response of Chamber Temperature to Reactor Power 

The frequency response of chamber temperature to reactor power at four steady - 
state operating levels is shown in figure 4. Data for these curves were obtained from the 
analog computer simulation with the turbine -control -valve position held constant at each 
level. 

The transfer functions that approximate these frequency -response curves are listed 
in table n . No attempt was made to derive these transfer functions from the system 
equations. The steady-state gain terms indicate an inverse relation to core weight flow. 
Attempts to correlate the time constants with any of the system variables were unsuccess- 
ful. 


TABLE H. - TRANSFER FUNCTIONS OF CHAMBER 


TEMPERATURE TO REACTOR POWER 


Steady -state 
operating 
level 

Transfer function, jATg(s)y^\Q(s)J 

9 V 

I 


2.6X10" 3 (l + ) 

V 1.55/ 


l + l- 4 s + /_s\ 2 
0.5 \0. 5/ 


n 


3. 67xl0' 3 (l + s ) 
V 1.55/ 


1 + 1- 4 s + / 8 ) 2 

0.45 \0. 45/ 

( 1+ m) 

m 


6. 75X10 -3 (l + s ) 
\ 1.25/ 


L 1 - 48 ^ s f 

0.27 \0. 27/ 

H) 

TV 

6.0X10" 3 ^1 + 


i* 1 - 48 */ s -f] 
0.26 \0.26/ 

(■^) 


Response of Chamber Pressure to 

Turbine-Control-Valve Position 

The frequency response of chamber 
pressure to turbine -control-valve posi- 
tion at the four steady -state operating 
levels is shown in figure 5. These data 
were obtained with constant control-drum 
position at each level. Approximating 
transfer functions for these frequency - 
response curves are listed in table III. 

Response of Reactor Power to 
Control -Drum Reactivity 

The frequency response of reactor 
power to control-drum reactivity at the 
four steady -state operating levels is 
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Phase shift, deg Amplitude ratio. 



(b) Phase shift. 


Figure 5. - Frequency response of reactor -out let-chamber pressure to turbine- 
control-valve position, obtained from analog computer simulation. Control- 
drum position held constant at each steady-state operating level. 
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table m. - transfer functions of chamber shown in figure 6. Data for these curves 

pressure TO turbine -control-valve position were obtained with the turbine- control - 

valve position held constant at each level. 

Designing reactor control systems 
requires a reactor transfer function dif- 
ferent from the one normally obtained 
for steady-state power operation. Con- 
trol systems designed by considering the 
steady -state reactor transfer function 
have exhibited unexpected oscillations 
during rapid power increases. Investi- 
gation of this phenomenon led to the ex- 
planation of the "period equilibrium 
effect" in reactor dynamics and the deri- 
vation of a reactor kinetics transfer func- 
tion that is dependent on reactor period 
(ref. 8). 

The data for the frequency -response 
curves shown in figure 6 were obtained 
by perturbing the system at steady-state 
power levels. Thus, the curves repre- 
sent essentially the response at infinite 
reactor period. The frequency response at finite -equilibrium periods is not readily ob- 
tainable from the analog computer simulation. Thus, another approach was used to ob- 
tain period -dependent reactor transfer functions. This approach is based upon the rela- 
tion among the reactor kinetics, the inherent reactivity feedback, and the desired reactor 
transfer function AQ/ASK^, as shown in the block diagram in figure 7. 

In appendix C, a transfer function that represents period-dependent reactor kinetics 
is discussed. Also presented in this appendix are specific transfer functions calculated 
for the system design values of l* and /3, and for infinite, 2- and 0. 5-second periods. 

The transfer functions of inherent reactivity feedback to reactor power for the four 
steady-state operating levels are listed in table IV. These transfer functions were de- 
termined by approximating the frequency response obtained from the analog computer 
simulation. 

The transfer functions listed in table IV were combined with the corresponding 
reactor -kinetics transfer functions for each operating level to obtain the desired closed- 
loop transfer functions AQ/ASK^. The ESIAC computer was used to find the roots of the 
equations that represent the denominators of the closed-loop transfer functions. The 
closed-loop transfer functions thus obtained are listed in table V. 
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Phase shift, deg Amplitude ratio, -j=±L 



Figure 6. - Frequency response of reactor power to control-drum reactivity, ottained from analog computer 
simulation. Turbine-control-valve position held constant at each steady-state operating level. 



Figure 7. - Block diagram of reactivity loop. 





TABLE IV. - TRANSFER FUNCTIONS OF INHERENT 


REACTIVITY FEEDBACK TO REACTOR POWER 


Steady -state 
operating 
level 

Transfer function, [a6 Kj(s)y^Q(s)J 

9 v 

I 

1. 4xl0~ 8 

(l + * )( 1+ Jl\ 
V 0. 48/ \ 71/ 

n 

2.1X10" 8 

(**£)( *5/ 

3.4X10 -8 

m 


IV 

3.27X10" 8 

(l + — )(l +— ) 
V 0.2/V 90/ 


The effect of reactor period on the reactor-power dynamic response is shown in fig- 
ure 8. In this figure are the frequency-response curves that represent the three period- 
dependent transfer functions for operating level I listed in table V. Also plotted in fig- 
ure 8 are frequency -response data obtained from the analog computer simulation at oper- 
ating level I. The agreement between analog data and the infinite period curve is reason- 
ably good. 
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TABLE V . - TRANSFER FUNCTIONS OF REACTOR POWER TO CONTROL -DRUM REACTIVITY 


Steady- Reactor Transfer function, |aQ(s)/4 5 K D (s)J 

state period, 0 v 

operating sec 
level 



Steady- Reactor 
state period, 

operating sec 
level 

m Infinite 












CONCLUDING REMARKS 


A nonlinear, wide-range, analog computer simulation was programed from a set of 
equations that represents the principal dynamics of a nuclear rocket engine system. 

From this simulation, the following frequency responses and transfer functions, for four 
selected steady-state operating levels, were obtained: chamber temperature to reactor 
power, chamber pressure to turbine-control-valve position, and reactor power to 
control-drum reactivity. 

The effect of reactor period on the transfer function of reactor power to control -drum 
reactivity was investigated. Transfer functions were determined for constant periods of 
2.0 and 0. 5 second as well as for an infinite period. 

The system transfer functions that were developed in this study have been used to 
perform a controls analysis (ref. 10) and to determine controller transfer functions and 
closed-loop responses of the system variables. 

In reference 10, the control system determined by root-locus analysis was imple- 
mented on the analog computer and used to control the system simulation. Closed- loop 
responses obtained from the analog simulation substantiated the analytically predicted 
responses. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 22, 1966, 

122-29-03-06-22. 
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APPENDIX A 


A area, ft^ 

11* 

C i concentration of 1 group of de- 

layed neutrons (see eq. (7)) 

c specific heat, Btu/ (lb) (°R) 

c specific heat at constant pressure, 

P Btu/ (lb) (°R) 

D diameter or hydraulic diameter, ft 

Fj( ) function of ( ), where i is an 
integer 

f friction factor 

2 

g gravitational acceleration, ft/sec 

H enthalpy, Btu/lb 

h convective heat-transfer coeffi- 

cient, Btu/ (ft^) (sec) (°R) 

I turbopump moment of inertia, 

(ft) (lb) (sec 2 ) 

J mechanical heat equivalent, 

ft-lb/Btu 

Kj constant, where i is an integer 

k thermal conductivity, 

Btu/ (ft) (sec) (°R) 

L length, ft 

l * mean effective lifetime of neu- 

trons, sec 

M torque, (ft) (lb) 

N turbopump speed, rpm 

n relative neutron level 

P total pressure, psia 


SYMBOLS 

p static pressure, psia 

Prandtl number, cu/k 

lr 

Q generated thermal power (o=n), 

Btu/sec 

q heat stored, Btu 

R gas constant for hydrogen, ft/°R 

Re Reynolds number, WD/Ap. 

s Laplace transform variable, sec' 

T temperature, °R 

t time, sec 

u velocity, ft/sec 

V volume, ft^ 

Q 

v specific volume, ft /lb 

W weight flow rate, lb/sec 

Cj divided by n 
/3 total delayed neutrons 

tVi 

/3 i fraction of i n group of delayed 

neutrons 

y ratio of specific heats 

A change or perturbation 

6K total reactivity 

6Kj inherent reactivity feedback 

6 angular position, deg 

X. decay constant of i n group of de- 

1 _i 

layed neutrons, sec 
[i absolute viscosity, lb/(ft)(sec) 

O 

p density of propellant, lb/ft 


Pr 
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t reactor period, sec 
Subscripts: 

A first reactor-core section 

a ambient 

B second reactor-core section 

b bulk conditions 
C third reactor- core section 

c reactor- core coolant passages 

cb cold bleed 
D control drum 
d demand 
e error 
H hydrogen 
hb hot bleed 
hg hot gas 
in inlet 
m material 
me reactor-core material 
mn nozzle material 
mr reflector material 
nc nozzle coolant passages 
nt nozzle throat 
o constant or steady-state value 
out outlet 
p pump 


r reflector coolant passages 

s surface 

t turbine 

ti turbine inlet 

te turbine exit 

tn turbine nozzle 

v turbine control valve 

w wall 

1 hydrogen tank 

2 pump inlet 

3 pump outlet 

4 nozzle coolant-tube inlet 

5 reflector inlet 

6 reflector outlet 

7 reactor- core inlet 

7. 1 first core-section outlet 

7. 2 second core- section outlet 

8 reactor outlet chamber 

9 bleed mixing volume 

10 turbine- control- valve inlet 

11 turbine inlet plenum 

12 turbine outlet plenum 
Superscripts: 

(~) average 

(') measured variable 



APPENDIX B 


SYSTEM EQUATIONS 

The specific system equations in this appendix are based on design data for a first- 
generation NERVA-type nuclear rocket engine that develops approximately 56 000 pounds 
of thrust. 


Pump 

P 2 = Pj = Const 
T 2 = Tj = Const 



T 3 = T 2 + °. 14 W p 
H 3 = 5. 0 T 3 - 298 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 


Pump Discharge Line 


T 4 = T 3 

(B7) 

h 4 = h 3 

(B8) 

dW r _ 4 . 2l 

p = 31. 57 P 3 - P 4 - 6. 56x10 (W^ 

(B9) 
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Nozzle 


p„ = p k + w: 


! / 8. 63xl0 -<s , 1. 81xl0“*\ 

(BIO) 

' \ p 4 P* / 


P 4 = Const 

(BH) 

P* 


. = 0. 1878 — - 

(B12) 


d<lmn = O. 32 8 (T 8 - T mn ) - 1. 035 W°; 8 ( T r ^ + ^ 


dt 


nc \ mn 


= F s(cu,«) 

mn o ^mn 7 


T c = 0. 294 

D 


H 5 + 147 - 1. 285X10' 5 (400 - P 5 )(240 - Tg) 


(B13) 

(B14) 

(B15) 


H 5 = H 4 + 1. 035 W-°- 2 ( Tmn - 


(B16) 



Reflector 


P 5 = P 6 + *?I 

fl. OlXlO -3 | 4. 12xl0" 3 \ 

(B17) 

\ p 5 p 6 / 


W„ = W n 
r p 

(B18) 

dVr = K^Q - 
dt ^ 

h r A s r(^mr - Ts + T0 ) 
r s,ry mr 2 I 

(B19) 

= F 4 (q) 
mr 4^mr 

(B20) 
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H 6 = 3.4T 6 - 


dq 

4 mc 

dt 


T e = 2 


W_ 


mr 


h r^s, r 


(h 6 - h 5 ) 


- T ( 


A = 87. 11 8 (T rnr )"°‘ 55 F, 

r s, r r ' mr' o 




Hg = 3. 85 Tg - 255 for Tg ^ 240° R 


(B21) 

(B22) 

(B23) 


147 + 1. 285x 10” 5 (400 - P g )(240 - Tg) for Tg < 240° R (B24) 


Reactor-Core Inlet Plenum 

dp^ • • • 

•^=0.4(W r -W c -W cb ) 

Pg = 5. 322 PgTg 

P 7 = P 6 
T 7 = T 6 

H 7 = H 6 

P7 = Pq 


(B25) 

(B26) 

(B27) 

(B28) 

(B29) 

(B30) 


Reactor-Core Section A 

- ' K A« - h c,A A s,c,A » 31 > 

Va'^Va’ (® 32 > 
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T 7 . 1 = 2 


W. 


me, A , * 

L n c, l A s, c, . 


(H? t - h 7 ) 


- T„ 


h c, A A s,c, A = 72 * 0 W®- ^mc, A>'°* 55f 7 ( T? 4 a T? - * ) 


H ? j = 3.51 T ? j - 90 


Reactor-Core Section B 


dq 


me, B 


dt 


k b Q " h c, B A s, c, B ( T mc, B 


T 7 . 1 + T 7 . 2 \ 


7 mc,B ^ 6 ^mc, 


T 7.2 " 2 


W. 


me, B . A 

n c, B a s, c, B 


^ 7.2 " H 7 . P 


- T 


7 . 1 


h s, B A s, c, B = 72 ‘ 0 V mc , B )"°- 55 F 7 p^^ 


Hr, 0 = 3.51 Tr, o - 90 for T„ „ < 2000 u R 


* 7.2 


■ 7.2 


7.2 


H ? 2 = 3 . 77 T 7 2 - 610 for 2000 ° R < T ? 2 ^ 3000 ° R 


H ? 2 = 4 . 09 T 7 2 - 1570 for T ? 2 ^ 3000 u R 


Reactor-Core Section C 


dq, 


" K cQ ■ h c, C A s, c, C ( T mc, C 


T 7. 2 + T 8 \ 


(B 33 ) 

(B 34 ) 

(B 35 ) 

(B 36 ) 

(B 37 ) 

(B 38 ) 

(B 39 ) 

(B 40 ) 

(B 41 ) 

(B 42 ) 

(B 43 ) 
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T mc, C = F 6^mc, Cp 


Tg = 2 


W, 


me, C h A 

n c, C A s, c, C 


(h 8 - h 7-2 ) 


- T, 


7.2 


h c, cv c, c - ,2 - ' ° v mc> c>'°‘ 


Ho = 3. 51 To - 90 for T 0 =s 2000 u R 


Hg = 3. 77 Tg - 610 for 2000° R< T g < 3000° R 


Hq = 4. 09 T q - 1570 for T c S 3000 R 


Reactor-Core Weight Flow 


w = 

C 


P 7 ' P 8 


v0. 5 


^4.475xl0" 3 v 7 + 3.807 x10‘ 4 v 8 


v 7 = 


P 7 


Vo = 


*>8 


Thrust Chamber and Nozzle 


dt 


= 0. 1(W C - W nt - W hb ) 


p g - 5. 322 PgTg 


(B44) 

(B45) 

(B46) 

(B47) 

(B48) 

(B49) 

(B50) 

(B51) 

(B52) 

(B53) 

(B54) 
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p ft 

W , = 8. 139 B 
nt T 0. 5 

o 

(B55) 

Reactor Kinetics 


6 

= (40 0006K - 312. 0)n + ^ \.C. 

i=l 

(B56) 

dC n 

1 = 11.856 n - 0.0127 C , 
dt 1 

(B57) 

dC 9 

* = 66.456 n - 0.0317 C 0 
dt 1 

(B58) 

dCo 

6 = 58.656 n - 0. 115 C, 
dt 6 

(B59) 

dC 4 

% - 126.984 n - 0.311 C, 
dt * 

(B60) 

dCe 

dt 5 

(B61) 

dC R 

dt 6 

(B62) 

Reactivity Summation 


6K = 6K d + 5K 0 

(B63) 

6K I = F 8^mc ,p H,c^ 

(B64) 
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p H,i 


P 7 T fi 

= o. 1878 1 — In — 

T 8 " T 7 T 7 


(B65) 


T * + T rj + T p 

m __ me, *x me, B me, C 

A mc o 


(B66) 


Bleed System 


dpj 

dt 


= 1.0(W cb+ W hb -W v ) 


Pg - 5, 322 PgTg 


dH 9 _ W hb H 8 , W cb H 6 W v H 9 


dt 


p 9 


P 9 p 9 


(B67) 

(B68) 

(B69) 


T 9 = 


Hg + 90 
3.51 


(B70) 


W b = 0. 542 AcbI>6 Fg( — 

t2- 5 VW 


(B71) 



(B72) 


(B73) 


\ ~ F 10^v^ 


P 10 ~ P 9 


(B74) 

(B75) 
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T 10 - T 9 


Turbine 


dp 


11 


dt 


= 1.0(W V -W ti ) 


P 11 " 5 - 322 P 11 T 11 


T 11 " T 10 


w - P 11 r /M2 1 
0 5 p 

Tjj V 11/ 


r 


M t = wJ 12 - 61 T ?i 5 


1 - 


' p 12 


0.2861°- 5 


11 / 


- 2. 85x10” J N 


Tj2 ~ T — ®35xl0 


_ 5 m^n 

w ti 


dp 


12 


dt 


= i. °(w ti - w tn ) 


P 12 - 5. 322 P 12 T 12 


= P i0 - 1.324X10 


p 12 “ *12 


_ w 2 

-2 _ti 

p 12 


A. P 10 /P o 

W tn = 0. 542 tn 12 F 0 [ — a 

T°2 5 \ 12> 


(B76) 

(B77) 

(B78) 

(B79) 

(B80) 

(B81) 

(B82) 

(B83) 

(B84) 

(B85) 

(B86) 
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Torque Balance 


— = 79. OOVL - M J 
dt F 




(B87) 


34 



APPENDIX C 


REACTOR-KINETICS TRANSFER FUNCTION 

The reactor-kinetics equations can be written in logarithmic form: 


— (In n) 
dt 


6K - /3 


Ev, 


dY. 

i 


dt 


~ t - x i y i 

r 


- Y A (In n) 
1 dt 


(Cl) 


(C2) 


The following transfer function can be derived (ref. 8) from these equations by using a 
perturbation analysis. 


A In n(s) _ 
A5K(s) 


l* s 


1 + 


£ 




i l (s + 


(C3) 


In the perturbation analysis that led to this transfer function, it was assumed that the 
steady- state power level n Q was constant and that the steady- state reactivity 6 K q was 
zero. 

Reference 8 shows that another reactor transfer function can be derived by assuming 
that the reactor power level is rising exponentially and that the reactivity is constant. 
This transfer function is 


A In n(s) 
A6 K(s) 



(C4) 


In reference 9, the following values of /T//3 and X i for a three- delay- group approx- 
imation of the reactor-kinetic parameters are given: 
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Xj = 0. 0276 


0! 

— = 0.274 

(3 


@2 

— = 0. 559 

0 


0Q 

- 0. 167 

0 


— 6 

By using these values, Z* = 25x10 second, 
tions were calculated for constant periods o 


X 2 = 0.235 
Ag = 1. 65 

and /3 = 0.0078, the following transfer func- 
infinity, 2 seconds, and 0. 5 second. 


For t q = °°, 


10.32(1 +- — - — ) (l + — 1 +- 
A In n(s) _ \ 0. 0276/ \ 0. 235/ \ 1.65, 

/ \ / \ / 

s 


A6K(s) 


s ( 1 + 


0. 095 


1 + 


1. 398 


1 + 


312.42 


(C5) 


For t q = 2 seconds, 


3. 848X10 2 (l + — \ (l + — \ (l + —^-4 
A In n(s) = \ 0. 528/ \ 0. 735/ \ 2. 15/ 

A6K(s) s L + L + _s__\ 4 + s ) 

\ 0. 541/ \ 1. 578/ \ 101.52/ 


(C6) 


For t q = 0. 5 second, 

2. 465xl0 3 ( 1 + — — \ (l + — ^ (l +— — \ 

A In n(s) V 2. 028/ \ 2. 235/ \ 3.65/ ( C 

A6K ( S > _ b (i + -L)(i + -*2 \(i + _JL_V 

\ 2.04/V 2.86/ \ 46.1/ 

The equivalent linear form of these transfer functions (i. e. , AQ(s)/A6K(s), with 
power expressed as thermal power) can be obtained by multiplying the logarithmic form 
by Q q , which is the steady -state power level. 
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